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We study the fluctuation-dissipation theorem for a Brownian particle driven into a nonequilib- 
rium steady state experimentally. We validate two different theoretical variants of a generalized 
fluctuation-dissipation theorem. Furthermore, we demonstrate that the choice of observables cru- 
cially affects the accuracy of determining the nonequilibrium response from steady state nonequi- 
librium fluctuations. 
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According to Onsager, in thermal equilibrium the re- 
action of a system to a small external perturbation and 
the decay of an internal fluctuation created by thermal 
noise are indistinguishable. This property, characteristic 
for the linear response around equilibrium, is expressed 
by the fluctuation-dissipation theorem (FDT) 



(1) 



which relates the time-dependent response Ra,h{t) of an 
observable a{t) to a perturbation h to the correlation 
function between a{t) and the derivative of the energy 
rate E with respect to /i Here, fceT is the thermal 
energy. 

Since the FDT allows to determine response properties 
like mobilities or susceptibilities from equilibrium mea- 
surements of diffusivities or power spectra and vice versa, 
it has found widespread application in different scientific 
fields like statistical mechanics, biophysics, chemical or 
solid state physics (H . In its original derivation the FDT 
holds only close to thermal equilibrium. Therefore, vio- 
lations of the FDT are a clear fingerprint of a nonequilib- 
rium system. So far, recent theoretical and experi- 
mental progress @ demonstrated that the FDT can also 
be extended to a specific class of nonequilibrium systems, 
i.e., nonequilibrium steady states (NESSs). In [8| it has 
been shown that a similar expression as in Eq. ([1]) can 
be obtained when the energy is replaced by the entropy 
within the correlation function on the right hand side. 
Furthermore, it was demonstrated that additional equiv- 
alent forms of the FDT exist, which in principle allows 
for infinitely many variants. 

In this Brief Report, we experimentally demonstrate 
the validity of the FDT in a NESS for two different vari- 
ants. Although both are equivalent from a theoretical 
point of view, large differences regarding the size of ex- 
perimental errors exist. Therefore, the right choice of 
observables is important for the accurate determination 
of the response in such measurements. 

Our experimental setup has been described in detail 
elsewhere ^^^\ and will be discussed here only briefly. 




FIG. 1: (Color online) Normalized steady state probability 
distributions Pu{x) (unperturbed state; solid red line) and 
Pp{x) (perturbed state; dashed black line). Inset: (a) Cor- 
responding tilted potentials of unperturbed (red solid line) 
and perturbed (black dashed line) state, (b) {sm{x/R)) for 
different driving forces. The black dashed line is a linear fit. 



A colloidal silica particle of radius r = 0.65 |J.m, im- 
mersed in water, is trapped within a three-dimensional 
torus of radius R = 1.18 |J.m by means of scanning laser 
tweezers. Using video microscopy, we track the angular 
coordinate x of the particle with —ttR ^ x < ttR with 
a spatial and temporal accuracy of 10 nm and 15 ms, 
respectively. Since the torus is far away from the lower 
surface of the cuvette cell (approximately 50 M-m), hy- 
drodynamic interactions with the sample cell are negli- 
gible [lij . The tweezers scanning motion exerts a non- 
conservative force / to the particle so that it circulates 
around the torus onto which an additional static sinu- 
soidal potential V{x) = (Vo/2) sin(a::/i?) is imposed. The 
corresponding total force 

F{x) = -d,V{x) + /, (2) 

with Vo = 58 keT and / = 30 kBT/^im [III, drives the 
particle into a NESS with an average circulation velocity 
of Tj*' = 7 [xm/s. In this NESS, the particle permanently 
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dissipates heat into the surrounding media, while the sys- 
tem is stiU characterized by a time-independent proba- 
bihty distribution p^{x), measured as shown in Fig. [U 
and a constant current j^. Brownian particles driven by 
forces such as Eq. (2) have been discussed in a variety of 
fields as reviewed in [TBt - 

As in thermal equilibrium, the response of a NESS 
is defined as the functional derivative of {a{t)) with re- 
spect to the perturbation h in the limit of h — > 0: 
Ra,h{t) = ^sh(Q) ■ However, in experiments such a func- 
tional derivative is not accessible. For a step-like pertur- 
bation of height A/i integration of Ra,h{t) leads to the 
integrated response 

Ra,h\T) dr = 



-Ah 



(3) 



which is experimentally accessible. Here, the average 
{o-(t))u is determined after the perturbation is turned 
off and {a)p is the average of a in the perturbed state, 
determined via a separate stationary measurement. In 
principle, the FDT allows for many choices of the ob- 
servable a and the perturbation h. However, due to 
the symmetry of our system, a convenient choice for the 
observable is a = sin(x/i?), which is obtained immedi- 
ately by the measured angular coordinate x. For the 
perturbation we choose Ah = A/, i.e., a small vari- 
ation of the driving force /. In a dynamical experi- 
ment, we switch between the reference NESS and the 
perturbed NESS with a period of 6 s for approximately 
4000 times by instantaneously reducing the driving force 
/ by A/ = —4 keT/iJ-m, while all other system param- 
eters are kept constant. The linear dependence between 
a and /, as shown in Fig. [Ijb), ensures that even for a 
force perturbation of 15 % the system stays within the 
linear response regime and that A/ = —4 keT/iJ-m is 
sufficiently small. 

First, we want to demonstrate the violation of the equi- 
librium FDT as given by Eq. ([T]). Since x = OajE is 
the conjugate observable to A/ with respect to E, we 
get ksT Ra.Af{t) = {o-{t)i{0)) with x the actual particle 
velocity. Thus, by defining for any observable b the in- 
tegrated correlation function as Cb{t) = (a(T)b(0))dr, 
the integrated equillibrium FDT following from Eq. ([T]) 
reads 

kBTnit)=Ci.{t). (4) 

Figure (Sl^a) shows kBTTZ{t) (open black circles), from 
a measurement involving the step-like perturbation, and 
the integrated normalized correlation funcion (red solid 
line), from a separate stationary measurement. Both 
quantities oscillate with a period of 1 s corresponding 
to the mean revolution time of the particle. This oscil- 
latory behavior is an inherent feature of this NESS [l^l , 
still present even though the motion of the particle is 
overdampcd and inertia is completely negligible. In con- 
trast, the response of an overdampcd system at thermal 
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FIG. 2: (Color online) (a) Integrated response kBTTZ{t) (open 
black circles) and integrated normalized correlation function 
C±(t) = C±{t) - {a){x)t (red solid line), (b) ksTTlit) (open 
black circles) and integrated correlation function C±_ys(t) 
(black solid line) according to the integrated generalized FDT 
[Eq. (O]. (c) Conditional steady state probability distri- 
bution p(i|a::o) at the flattest a::o = —1.9 li-m (black closed 
bars; u^{xo) ~ 3 Virn/s, F{xo) — 8 kBT/jim) and the steep- 
est xo — —2.6 [im (red open bars; i''^{xo) = 20 (xni/s, 
F{xo) — 55 kBT/[j.m) position in the tilted potential. 



equilibrium always decays exponentially. The fact that 
the response is more than one order of magnitude smaller 
than the correlation clearly demonstrates that the system 
is far away from thermal equilibrium. 

The FDT can be generalized to a NESS by choosing 
the observable as conjugate to the perturbation A/ with 
respect to the system entropy s{t) = — fee \iip'^{x{t)) Q, 



kB RaAfit) = (a(t)h9A/s(0)]). 



(5) 



In thermal equilibrium Eq. ([5]) becomes Eq. ([T]) since in 
the absence of an external driving force 9a/ s reduces to 
dAfE/T It is convenient to split the entropy pro- 
duction of the system into one of the medium Smed and 
a total entropy production stot [l^ . 



Stot 



(6) 



According to Clausius, s,ncd = Q/T ~ xF/T corresponds 
to the heat flow Q into the thermal bath. The total 
entropy production rate is given by stot = xv^/{^qT) 
with the bare mobility and the local mean velocity 
i'^[x) = j^/p^{x), which corresponds to the average over 
all stochastic velocities for a fixed particle position x. 
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After inserting Eq. ^ into Eq. ((5)) and integrating over 
time, the latter becomes 0, Q 



0.05 



kBTTZ{t) =C±{t) -C^s{t) = C±-y.{t). 



(7) 



Here, the first correlation function corresponds to the 
one occurring also in the integrated equilibrium FDT 
[Eq. (HI] but now evaluated under nonequilibrium con- 
ditions. The second term reflects the violation of de- 
tailed balance, which vanishes under equilbrium condi- 
tions iy'' 0). 

Our measurements [see Fig. [^b)] confirm the general- 
ized FDT from Eq. ([7]) connecting response and correla- 
tion function since within the statistical errors kBTTZ{t) 
(open black circles) and C±-,^o (t) (black solid line) coin- 
cide. We want to emphasize that Eq. ([7]) was tested here 
over more than two relaxation times (r^ = 1.4 s) corre- 
sponding to about three particle revolutions. Due to this 
short relaxation time (in our experiment ^ R^/Dq 
with Do = k^T / (%'Kirjr) the free diffusion coefficient and 
■q the viscosity of the solvent [13]), in contrast to pre- 
vious work [9|, we were able to observe the generalized 
FDT along the full damped oscillations of kBTTZ{t). 

The FDT in the NESS can be brought into a form even 
more reminiscent to the equilibrium FDT by introducing 
the velocity v{t) = x{t) — i'^{x{t)), which is measured 
with respect to the local mean velocity 0, 0, 0] . Equa- 
tion d?]) then reads k'QTTZ{t) = Cv{t). A physical reason 
for this restoration of the equilibrium FDT in the locally 
comoving frame can be seen in measuring the conditional 
velocity distributions p{x\xo) in this frame. Their shapes 
are the same at any position xq as illustrated in Fig.[2l[c) 
for two specific positions. Both distributions resemble 
Gaussians with a mean equals to the local mean velocity 
ly^i^xo) and a width of = 2Do/At, where At = 20 ms 
is the time interval over which the velocity x was deter- 
mined. Thus, locally the velocity fluctuations are indis- 
tinguishable from equilibrium ones even though nontriv- 
ial velocity correlations develop in a NESS even when 
expressed in the locally comoving frame. 

Somewhat surprisingly, the conjugate observable ap- 
pearing in the correlation function on the right hand side 
of the generalized FDT in Eq. ^ is not unique as ob- 
served in 0, where a classification of the different vari- 
ants was suggested. This equivalence in a theoretical 
perspective, however, does not extend to an equivalence 
in terms of experimental accessibility as we next demon- 
strate. The earliest variant of the FDT in a NESS is 
based on Agarwal [l3] and reads in our notation 



kBTRa.Afit) - {ait)[~Dod,lnp'{x{Om. 



(8) 



Here, the conjugate observable is derived from the steady 
state distribution which leads to an observable in the 
state (i.e., configuration) space that does not involve a 
time derivative. Consequently, it is more easily accessible 
experimentally than the strongly fluctuating stochastic 
velocity x appearing in Eq. ([T])- 




FIG. 3: (Color online) (a) Integrated correlation functions 
Ci_^s(/:) (black dashed line) and Ci/=-^i(,F(t) (red solid line) 
indicating different forms of the generalized FDT. (b) Com- 
parison between kBTTZit) (black open circles) and Cys-^QF(t) 
(red solid line). 



In a steady state, the constant probability current can 
be written as f = fioF{x)p''{x) - Dodj^p^x) With 
the definition of the local mean velocity, we obtain 

- Dod,: lnp%x) = ly'ix) - ^^oFix). (9) 

Inserted into Eq. (|8]) and integrated over time, we get 

kBTUit) ^ C,.it) - C^„F{t) = C^s^^oFit). (10) 

Since both the integrated generalized FDT from Eqs. ([7]) 
and pO)) are valid for any steady state, their correlation 
functions must be identical 
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Even though their second arguments are different, there 
is indeed a very good agreement between these correla- 
tion functions, as seen in Fig. EJa). The deviation, which 
occurs only for times larger than 1.5 s, can be attributed 
to statistical errors which increase in time. Compared to 
the previously shown function C±-vs{t) [see Fig.[2][b)] the 
correlation function C^=_^j,i?(t) [see Fig. |3{b)] traces the 
integrated response even better. As most pronounced 
in the time interval between 1.5 and 3 s, the phase 
and the amplitude of the oscillations are represented by 
Ciy=-pf,F(i) more precisely. 

To understand why C,jB_^gp{t) leads to a more accu- 
rate determination of kBTTZ{t) we next concentrate on 
the terms involved in these two different correlation func- 
tions. Figures SJ^a) and SJ^b) show their time evolution 
along the trajectories over a time interval of 3 s. Al- 
though leading to the same correlation function when 
correlated with a = sin(a;/i?) (light gray line), the tra- 
jectories of — hqF (red line) and x — v^ (black line) are 
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FIG. 4: (Color online) (a), (b) Quantities involved in the cor- 
relation functions: sin(a::/i?) (light gray line), — fj,oF (red 
line) and x — (black line), (c) Integrated correlation func- 
tion C±-i^a (t) determined from different raw data lengths. The 
black line corresponds to the full data set of 3000 s whereas 
for the blue dashed (red dashed-dotted) line only one half 
(quarter) of the original data was used. 

different. While the first trajectory is rather smooth, in 
the second one, due to the appearance of i, strong fluc- 
tuations are present. The amplitude of these fluctuations 
is more than ten times larger than the maximal velocity 
variation of — fJ-oF- Since for a reliable average (. . .) 
many reaUzations have to be sampled, it is evident that 
the correlation Cx-u" (t) cannot be determined with the 
same statistical precision than Cvo-^apit)- This observa- 
tion is supported by Fig. St^c) where, in order to examine 
the impact of statistics, only parts of the acquired raw 
data were evaluated. The black line shows Cx-u' (t) calcu- 
lated from the complete raw trajectory of length 3000 s. 
For the blue and the red line the length of the trajectory 
was reduced by a factor of two and four, respectively. In 
the latter case the strong deviation to the black curve is 
obvious. In contrast, C^n-fj^gpit) is almost insensitive to 
a reduction of the statistics. Indeed, reducing the length 
of the evaluated data to 300 s (one order of magnitude) 
does not influence the shape of C,^B_^gp(t) visibly. 

The insight gained in this model system can be trans- 
ferred to other systems as well: If one wants to obtain the 
response of a NESS from a measurement of its stationary 
correlations, one should use the variant of the generalized 
FDT not involving time derivatives for its better statisti- 



cal convergence properties. In this respect, the situation 
in experiments is somewhat different from that in simula- 
tions, where the correlation between a{t) and the a priori 
known Langcvin noise {i — fJ.nF)/2 = determines the 
response most conveniently [8|. 

In summary, we studied the FDT for a driven Brow- 
nian particle experimentally. For this paradigmatic sys- 
tem, a kind of " Ising-Model" of nonequilibrium steady 
states, we validated two different generalizations of the 
FDT to the nonequilibrium regime. We demonstrated 
that the right choice of observables affects the errors 
when calculating the response from a measurement of 
stationary correlations. 
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